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Abstract 



> 

\^ ', In this letter, we reevaluated the Higgs-mediated contribution to /x — )• 67, 

cn ' /U — >• 3e, and /i-e conversion in nuclei in the MSSM, assuming left-handed slep- 

l/^ ■ tons have flavor-mixing mass terms. Contrary to previous works, it is found that 

^P , Barr~Zee diagrams including top quark give dominant contribution to /x — )• 67, 

and those including bottom quark and tau lepton are also non-negligible only when 
tan/3 is large. As a result, the Higgs-mediated contribution dominates over the 
gaugino-mediated contribution at one-loop level in /i — )• 67 when Msusy/niAo ~ 50, 
rS ' irrespectively of tan/3 as far as tan/3 is not large. Here, MsusY and rriAO are a 

C^ ■ typical mass scale of the SUSY particles and the CP-odd Higgs boson mass, respec- 

tively. Ratio of branching ratios for // — )• 67 and /i-e conversion in nuclei is also 
evaluated by including both the gaugino- and Higgs-mediated contributions to the 
processes. It is found that the ratio is sensitive to tan/3 and MsusY/mAO when 
-^SUSyM^o ~ (10 - 50) and tan/3> 10. 



1 Introduction 

Charged lepton-flavor violating (cLFV) processes, such as /i — )■ 67, are sensitive to physics 
beyond the standard model (SM) [Ij. While the lepton-flavor conservation is not exact in 
nature due to finite neutrino masses, cLFV processes are quite suppressed in the standard 
model. Thus, searches for cLFV processes are a unique window to physics beyond the 
SM, especially, at TeV scale. 

Now the MEG experiment is searching for /i — )■ 67 [2], it would reach to ~ 10^^^ for 
the branching ratio on the first stage, which is improvement of two orders of magnitude 
compared with the current bound. The COMET and Mu2e experiments [3l H], which 
are searches for /i-e conversion in nuclei, are being planed in J-PARC and Fermilab, 
respectively. It is argued that they would reach to ~ 10~^^ for branching ratio of fi-e 
conversion with target Al. Here, branching ratio of /i-e conversion is ratio oi fi-e conversion 
rate over muon capture rate. Searches for fi ^ e'j and fi-e conversion in nuclei are 
complementary to each other in studies of physics beyond the SM since those processes 
may be induced by different types of processes. 

The minimal supersymmetric (SUSY) standard model (MSSM) is a leading candidate 
for physics beyond the SM, and cLFV processes are extensively studied in the model. 
SUSY-breaking slepton mass terms are lepton-fiavor violating. It is noticeable that ratios 
of branching ratios for cLFV processes would give information of mass spectrum in the 
MSSM, since dominant diagrams in cLFV processes depend on the mass spectrum. 

When SUSY particle masses are ^O(l) TeV, the muon LEV processes, such as /i-e 
transition processes, fi — )■ 67, yU — )■ 3e and /i-e conversion in nuclei, are generated by the 
gaugino-mediated contribution, which is generated by one-loop diagrams of gauginos and 
sleptons (and Higgsinos). Branching ratios for the cLEV processes due to the gaugino- 
mediated contribution are suppressed by l/M^^^y, since the effective dipole interaction 
is dominant in the cLEV processes. Here, Msusy is a typical mass scale of the SUSY 
particles. On the other hand, when Msusy ~ 0{l) TeV, the Higgs-mediated contribution 
to the processes could be sizable. The non-holomorphic LEV correction is generated to 
Yukawa coupling of the Higgs bosons at one-loop level, and it is not suppressed by Msusy 
[6]. Branching ratio of /i-e conversion in nuclei is more sensitive to the Higgs-mediated 
contribution [7]. Thus, ratio of branching ratios for /i — )■ e7 and /i-e conversion in nuclei 
is a good observable to constrain mass spectrum in the MSSM, since it is sensitive to 
whether the gaugino-mediated or Higgs-mediated contribution is dominant. 

In this letter we systematically calculate the Higgs-mediated contributions to cLEV 
reactions in the MSSM, and clarify the dominant process in each cLEV reaction. Eor 
this purpose, we first evaluate the Higgs-mediated contribution to /i — )■ e7 in the MSSM. 
Barr-Zee diagrams give dominant contribution to /i — )■ e7 among various diagrams though 
those are of higher order. We systematically evaluate those diagrams, and find that 
Barr-Zee diagrams including top quark give the largest contribution, and the branching 
ratio for /i ^ e7 induced by the Barr-Zee diagrams is approximately proportional to 
tan^ 13. The angle /3 is defined by tan/3 = (-^2) / (-^?)- The Higgs-mediated contribution 
dominates over the gaugino-mediated one in /i — )■ e7 when Msusy/'"^ao ~ 50, which is 



almost insensitive to tan/3. Here, m^o is the CP-odd Higgs boson mass in the MSSM. 

Using this result, we evaluate ratio of branching ratios for /i — )■ 67 and /i-e conversion 
in nuclei. When the Higgs- mediated contribution is dominant, the ratio of the branch- 
ing ratios is scaled by tan^/3. It is found that the ratio is quite sensitive to tan/3 and 
^susy/"^ao when Msusy/'^^ao ~ (10 — 50) and tan/3 ^ 10. We also check that the ratio 
of the branching ratios for /x — )• 67 and /i — )■ 3e is insensitive to them. 

In Refs. [SI [H] the Higgs-mediated contribution to the fi-e transition processes in the 
MSSM is discussed. It is argued that when the Higgs-mediated contribution is dominant, 
the Barr-Zee diagram including W boson is dominant and branching ratio of /i — )■ 67 is 
scaled by tan^/3, not tan^/3. This is obviously overestimated. We clarify what is wrong 
in their deviation. 

We assume that left-handed sleptons have flavor-mixing mass terms in this letter, 
simply because this setup is well-motivated from the SUSY seesaw model [5j. Extension 
to more general cases will be given elsewhere. 

This letter is organized as follows. In the next section we evaluate the Higgs-mediated 
contribution for /i — )• 67. We show ratio of the Higgs-mediated and gaugino-mediated 
contributions. In Sec. 3, we discuss the Higgs-mediated contributions to /i — )■ 3e and /i-e 
conversion in nuclei, and evaluate ratios among the cLFV processes. Sec. 4 is devoted to 
conclusions and discussion. 

2 Higgs-mediated contribution to /i — > 67 in the MSSM 

In the MSSM, LFV in the Higgs coupling originates from the non-holomorphic correction 
to Yukawa interaction of charged leptons |^j . By including the correction due to one- loop 
diagrams of gaugino and sleptons, the effective Yukawa coupling is given as follows: 

-£efr = e'^.ye^H',e',, + e'^.y,, (eP<5., + e^^^^) H^,*e',^ + h.c, (1) 

where yei stands for the i-th charged-lepton Yukawa coupling constant at tree level and 
e'^j and e^j represent right-handed and left-handed leptons, respectively, in a basis where 
the tree-level lepton Yukawa matrix is diagonal. The non-holomorphic interaction eg 
{i ^ j) is generated by flavor- violating slepton mass terms. As mentioned in introduction, 
we assumed that left-handed sleptons have flavor-violating mass terms. We parametrize 
€2 with mass insertions (Mis) parameters, 5^^^ = (Am? )ij/mj , where (Am? )ij is off- 

■^ I'L <'L ''L 

diagonal element of left-handed slepton mass matrix and m^- is an average left-handed 
slepton mass. When the SUSY-breaking mass parameters in the MSSM are taken to be 
a common value (Mgusy), the non-holomorphic corrections e^ and 63 are reduced to 
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Note that tf and e2 *■*■''' do not vanish even in a hniit of large masses of SUSY particles. 
This is quite different from LFV effective dipole operators induced by the gaugino-slepton 
loops, whose coefficients are suppressed by masses of internal SUSY particles. 

In a mass-eigenstate basis for both leptons and Higgs bosons, £eff for Z^-e transition is 
described as |6] 

- £^- = ^ ^" {jiPLe) [cos(« - /3)/i° + sin(« - (i)H^ - ^A°] + h.c. , (3) 

where h^ and H^ are the CP-even Higgs fields (rn^o < rriHo), and A^ is the CP-odd Higgs 
field. The LFV parameter A^g is given by A^g = £2/(1 + ^i tan/3)^. When we treat 
ef and 62 as a perturbation, we may neglect e-^'* of the denominator at the first order. 
In this letter, we set A^g = €2 ■ 

In the MSSM, LFV interaction of /i^ in Eq. ([3]) vanishes when the masses of H^ and 
A^ go to infinity, since cos(a — /3) behaves as 

cos(a-/3)~^^^. (4) 

This comes from a fact that SM does not have LFV and the light Higgs boson h'^ becomes 
SM-like in above limit. Therefore the contributions from H^ and A^ should be included 
in cLFV processes. 

Now we consider the Higgs-mediated contribution to /i — ?■ 67 in the MSSM. Effective 
amplitude for yU — )■ 67 is parametrized as 



T = ee*'^(g)lleb-g) 



m^ta^,q'{A'^PL + A''Pn 



u^^ip) , (5) 



and branching ratio of /i — )• 67 is derived as BR(/i — )■ 67) = {ASir^aem/ Gp){\A^\'^ + {A^l"^) . 
Here, aem(= e^/47r) is the fine structure constant and Gp is the Fermi constant. While 
this amplitude could be induced at one- loop level (Fig. [1]), it is suppressed by three chiral 
flips. I.e., one chirality flip in the lepton propagator and two lepton Yukawa couplings. 
Indeed two-loop diagrams may be significant contribution. As shown in Fig. [21 two-loop 
diagrams, called as Barr-Zee diagrams, involve only one chiral flip (from lepton Yukawa 
coupling), and hence their contribution is much larger than that at one- loop level. 

Following Refs. |S1 H] , we consider Barr-Zee diagrams which involve effective 7-7-0° 
vertices {(p^ = h^, H^, and A^). The effective vertices are induced by heavy fermion 
or weak gauge/Higgs boson loops. Barr-Zee diagrams involving bottom- and top-quark 
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Figure 1: /x — > 67 induced by Higgs boson exchange at one loop level. 
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Figure 2: Examples of two- loop Barr-Zee diagrams induced by Higgs exchange. 



loops (Fig. |2] (a)) give contributions to the coefficient A^ in Eq. ([5]) as 
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Here, Nc is color factor, and Qh(t) represents electric charge for bottom (top) quark, zlo = 
w? Jm^o for 0° = /i°, if°, A° and q = b, t. Similarly, we calculate the coefficient for tau- 
lepton loop by substituting A^^ = 1, and replacing Qbyrrih to Qr,mT-. The functions f{z) 
and g{z) are called Barr-Zee integrals, whose explicit forms and asymptotic behaviors are 
given in Appendix. For ?tt,^o ^ rri^o and tan/3 ^ 1, the Barr-Zee diagram contribution 



is approximated as 
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It is found that tan/3 and/or large logarithmic factors enhance heavy-Higgs (i/°, A^) 
contributions, and the light-Higgs {hP) contribution is subdominant. 

Similarly, the Barr-Zee contributions from loops of W~ boson (Fig. E] (b)), Nambu- 
Goldstone (NG) boson G~, and charged Higgs boson H~ are calculated. Each contribu- 
tion to A^ is derived as follows, 
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where z^o — '"-(/,-/ '"-^o 
from coupling of H^H'cfP , 
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/i°,i7°), and /^o (00 = /iO,/fO) comes 



ffiO = —2mly- sin(a — /3) + m|o sin (a + /3) cos 2/3 
/ji^o = 2m^_ cos(a — /3) — m|o cos (a + /3) cos 2/3 . 
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The functions F(z) and F'(z) are 3/(z) + 5c/(z) +3/4^(2) +3/4/1(2) and {g{z)- f{z))/{2z), 
respectively. h{z) is also given in Appendix. The CP-odd Higgs boson A^ does not 
appear here if CP is conserved. The above contributions are approximated in a limit of 
rriAO ^ rrizo and tan /3 ^ 1 as 



aR 

^BZ{W^-) 



ae.y2G^2m|o^^^^^, 



IGtt 



m 



AO 



/ie 



X 



F(C0S t^H/) — 



m 



AO 



log- 



m 



w- 



m 



AO 



(12) 



^BZ(G-) 



-acmV^Gi. 2m|o ^ .^ 
2^ tan /3 A 



167r3 



771 



AO 



L 
fie 






(13) 



^BZ(//-) 



ae. 



.V2G 



F"^|o 



167r3 



771 



tan/3A 



AO 



L 
fie 



X 



(Acos'ew^-2)F'{l] 



771 



(2 cos^ 9w-l 



771 



zo 



AO 



771 



AO 



log 



771 



AO 



i^o 



+ 



(14) 



Here, ^(cos^^vk) = 7.96, ^'(cos^^vk) = 0.121, and F'(l) = 0.172. 

Notice that the diagram of H*^ including W~ loop gives contribution to A^, which 
is proportional to tan/3 in a limit of m^o ^ tti^o and tan/3 ^ 1. This is because the 
H^W^W~ coupling is suppressed by cos(a — /3) ~ — 2m|o/(m^o tan/3) while the LFV 
Yukawa coupling in Eq. 03]) is proportional to tan^ /3. All contributions from Barr-Zee 
diagrams to Ar are proportional to tan/3 in a limit of m^o ^ m^o and tan/3 ^ 1. 
The exception is those including the bottom-quark and tau-lepton loops, which tend to 
be subdominant in moderate tan/3. It is argued in Refs. [H [9] that the diagram of H^ 
including W~ loop gives a contribution proportional to tan^ (3 and that it is the largest 
among the Barr-Zee diagrams. However, using corrected tan/3 dependence, this W~ loop 
contribution is no longer the largest one, as will be shown below. 

We also include Barr-Zee diagrams including H'^ loop, though the contribution is also 
subdominant. It is also proportional to tan/3 at most, and coefficient for its logarithmic 
term is suppressed by m\o/m\o- 
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Figure 3: Absolute values of coefficients for Higgs-induced dipole operator, |A^|, as a 
function of CP-odd Higgs boson mass 771^0- We show those from diagrams including 
heavy-fermion loops, diagrams of light- and heavy-Higgs bosons including W~ and G~ 
loops. For comparison, one-loop contribution to A^ is also shown. Here, left figure is for 
tan/3 = 50 and right one is for tan/3 = 10. We took A^^ = 5 x 10"*^, m^^o = 120 GeV, 
mt{mzo) = 181 GeV, and rrihirnzo) = 3.0 GeV. 



fie 



Fig. |3] shows each contribution to A^ as a function of CP-odd Higgs boson mass m^o 
for tan/3 = 50 (left) and tan/3 = 10 (right). Here, A^, = 5 x 10-^ m^^o = 120 GeV, 
mt{mzo) = 181 GeV, and mb{mzo) = 3.0 GeV. Barr-Zee diagrams including top-quark 
loop give dominant contribution to yU — )■ 67, and the bottom-quark one is also sizable only 
when tan/3 is large. The W~ and NG-boson diagrams tend to be subdominant unless m^o 



is small. It is noticed in Refs. [HI E] that there are strong cancellation between Barr-Zee 
diagrams of H^ involving W~ and G~ loops [TU]. However, other contributions dominate 
over them, and hence this cancellation effect does not appear in the branching ratio. 

For comparison, we show the contribution from the one-loop diagrams (Fig. [1]) in 
Fig. 131 It is approximated as 
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As expected, this is always subdominant. 

Now we consider competition between the gaugino- and Higgs-mediated contributions 
to /i — !■ 67. The gaugino-mediated contribution to A^ is approximated as 



" -i^S(i+i*-^«»-'W ^"' 



-^S.ugi.o = -n^;^(l + jtari^»„.)-rj5— ^S'™^. (16) 



where we take a common value Msusy for the SUSY particle masses. 

Fig. m is a contour plot for square of ratio of the Higgs- and gaugino-mediated con- 
tributions to A^ as functions of Msusy /""^ a" and tan/3. The line on which this ratio is 
equal to unity is boundary of the two regions where each effect dominates fi — )■ 67. In 
small tan/3 region, the Higgs-mediated contribution comes from mainly Barr-Zee dia- 
grams including top-quark loop. Both Higgs- and gaugino-mediated contribution to Aji 
are approximately proportional to tan/3. However, in large tan/3 region, the Barr-Zee dia- 
grams with bottom-quark and tau-lepton loops receive larger tan/3 enhancement, and the 
Higgs-mediated contribution becomes larger with the same Msusy /"^ao value. We found 
that the Higgs- and gaugino-mediated effects are comparable to each other in /x — )■ 67 
when Msusy/''^ao — 50. 

3 Correlation among LFV processes 

Now we discuss other /x-e transition processes, /x — )■ 3e and /i-e conversion in nuclei, when 
the Higgs-mediated contributions are dominant in the MSSM. These two processes have 
strong correlation with /i — )■ 67 when the gaugino-mediated contributions are dominant, 
since effective LFV dipole operator determines the processes. 

First, we consider /i — )■ 3e. This process is generated from three types of effective four- 
Fermi operators; scalar-, vector-, and dipole operators. Scalar operators are induced by 
tree-level Higgs boson exchange. On the other hand, vector and dipole ones are generated 
by virtual-photon mediating processes fi — )■ 67* at higher order. When only the Higgs 
bosons contribute to /i — )■ 3e, vector operator mainly comes from one-loop diagrams, and 
dipole one is generated by two-loop Barr-Zee diagrams. Since diagrams for vector and 
scalar operators need two chirality flips, these operators are suppressed by small Yukawa 
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Figure 4: Contour plot of square of ratio for Higgs- and gaugino-mediated contributions 
to A^ as functions of Msusv/mAo and tan/3. Higgs-mediated contribution is dominant at 
right-handed side of hnef on which this ratio is equal to unity. 
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couplings (y^ or t/e), compared to dipole operator. Vector and dipole operators come from 
higher-order effects and they are suppressed by loop factors. 

Thus, contributions of these operators are estimated roughly as follows, 
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where lower indices (0,1,2) mean coefficients for scalar, vector, and dipole operators, 
respectively. Ratio of these coefficients becomes Aq : Ai : A2 — I : 0(1) : 0(10), and the 
coefficient for dipole operator is the dominant contribution. There is also log {rn?^/rn?^ 
enhancement for dipole operator contribution to ;U — ;■ 3e, which comes from final state 
phase space integral. As a consequence, /x — )■ 3e is dominated by dipole operator and 
there is strong correlation between /i — )■ 3e and /i — )■ 67, 
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Next, we discuss /x-e conversion in nuclei. Dominant contribution to this process 
comes from tree-level Higgs boson exchange [7]. The reason is that coupling between 
Higgs boson and nucleon is characterized by the nucleon mass mjq through the conformal 
anomaly relation |rij , and could evade suppression of light constituent quark mass. Then, 



branching ratio for fi-e conversion in nuclei at large tan/3 is derived from formulae in 
Refs. Ha 113] as follows, 



BR(M1 ^ eAl) ~ 6.8 x 10^^ ^ ",," (A 
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Here, w^U - 



0.7054 X lO^sec"^, and nip is proton mass. We use the recent lattice 
simulation result [13] for the a term, which shows that the strange quark content of 
the nucleon is much smaller than previously thought. Notice that branching ratio for 
fi-e conversion in nuclei is scaled by tan^ (3, while those for fi ^ e'j and /i — )■ 3e are 
proportional to tan^/3. 
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Figure 5: Branching ratios of Higgs-mediated cLFV processes. We took tan/3 = 50 (left) 
and tan/3 = 10 (right), A^, = 5.0 x 10"^. 



In Fig. [5] branching ratios of Higgs-mediated LFV processes are shown as a function 
of TTi^o. Though we include contributions from the scalar and vector operators in the 
evaluation of BR(yU — > 3e) in addition to the dipole one, it is found from this figure that 
there is still tight correlation between /i — !■ 67 and /i — )■ 3e, BR(/i — )■ e7)/BR(yLi -^ 3e) ~ 
O(aem)- Thus, it is a signature that dipole operator dominates these two processes. On 
the other hand, fi-e conversion in nuclei is dominated by tree-level Higgs boson exchange, 
and such simple correlation does not appear, as expected. In the gaugino-mediation case, 
the dipole operator dominates three processes. Thus, it is important to measure /i-e 
conversion rate for discrimination of these two cases, in addition to /i — )■ 67. 

It is also found that while /i-e conversion process is simply scaled as 1/m'^o, other 
two processes are not. This is because other two processes receive large logarithmic 
corrections. 

Fig. [6] shows contour plot of BR(yuAl — > eAl) / BR(yU — )■ 67) including both the Higgs- 
and gaugino-mediated contributions. If the Higgs-mediated contribution is dominant in 
the cLFV processes, the ratio between /i — )■ 67 and fiN — )■ eA^ is sensitive to tan/3, but 
not to MsvsY/inAO- On the other hand, if gaugino-mediated LFV is dominant, this ratio 
is about 0{aem) since dipole operator contributions dominate the cLFV processes. When 
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Figure 6: Contour plot of BR(/iAl — )• eAl) / BR(yU — )• 67), tan/3 vs Msusy/''^ao including 
both the Higgs- and gaugino-mediated contributions. 



^susy/''tiao ~ (10 — 50) and tan /3;^ 10, both Higgs- and gaugino-mediated diagrams 
contribute to those processes in different way and we could give constraints Msusy/''^ao 
and tan/3 from BR(/xAl -)■ eAl)/BR(/i -)■ 67) . 



4 Conclusions and discussion 

In this letter, we reevaluated /x-e transition processes induced by non-holomorphic Yukawa 
interactions in the MSSM. We discussed correlation among branching ratios for /x — )■ 67, 
/i — )• 3e, and /z-e conversion in nuclei in the MSSM, by including both the gaugino- 
and Higgs-mediated contributions to the processes. It was assumed in this letter that 
left-handed sleptons have flavor-mixing mass terms. 

Though Higgs-mediated contribution to /i-e transition processes is evaluated in pre- 
vious works [HI [H], we found that contribution from Barr-Zee diagram including W~ 
boson, which was thought to be the largest contribution to yU — )■ 67 among various Higgs- 
mediated contributions, has incorrect dependence on tan/3. As a result, branching ratio 
for /i — )■ 67 was overestimated. We showed that Barr-Zee diagrams including top quark 
are rather dominant, and those including bottom quark and tau lepton are also sizable 
only when tan/? is large. Then, the Higgs-mediated contribution dominates over the 
gaugino-mediated one in /i — )• 67 when Msusy/rriAO ^ 50, irrespectively of tan/3 as far as 
tan/3 is not large. 

We evaluated ratio of branching ratios for /i — )■ 67 and /i-e conversion in nuclei by 
including both the gaugino- and Higgs-mediated contributions to the processes. We found 
that the ratio is sensitive to tan/3 and Mgusy/'^i^o when Msusv/mAO ~ (10 — 50) and 
tan/3 ^ 10. Ratio of the branching ratios for /i — )■ 67 and /i — )• 3e is insensitive to tan/3 
and the MSSM mass spectrum, since the dipole term contribution is always dominant in 
/i — !■ 3e. 
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In general, right-handed slepton mass terms or A terms could be sources for flavor- 
mixing. In particular, gaugino-mediated contributions from right-handed slepton mass re- 
ceive destructive interference between the bino and bino-Higgsino amplitudes ^5\. There- 
fore, Higgs-mediated contribution may be signiflcant in some parameter region and de- 
coupling behavior of Msusy could be modifled. We leave it for our future work. 

A Barr— Zee integrals 

The Barr-Zee integrals f{z), g{z) and h{z) are given by 
j{z) = -z ax — ^; ^log' 



(20) 



2 Jq x{1 — x) — z z 

1 f\ 1 x{l-x) 

h{z)(=z^^('-^ 
V dz\ z 



"^ dx 

2 Jq z - x{l - x) 



z , x(l — x) 

1 H r, V log ■ 



Z — x{l — X) z 

In the limit oi 1 ^ z, the asymptotic forms of them is given as follows [10 



f{z) ~ -{logzf , g{z) ~ -(log2;)^ , h{z) ~ zlogz , 
f{z)-g{z)r^zilogz + 2) . 



On the other hand, in the limit of 1 ^ z, 

1 13 1 1 

f{z) ~ -log 2; + —, g{z) ~ -log2; + l, h{z) ~ --(log2 + l) 

f{z)-g{z) r^--\ogz-—. 



35 3 

F{z) ~ —z{logzY + -zlogz 

F'{z)^-^{logz + 2) . 



On the other hand, in the limit of 1 <^ z. 



r./ N 7 181 

F{z) ^ -\ogz + — , 

F'(z)r~.— (- log z + —) . 
^ ' 2z^& ^ 1%' 
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(22) 



Similarly, in the limit of 1 ^ 2;, the asymptotic forms of F{z) = 3f{z) + 5g{z) 
3/4:g{z) + 3/4:h{z) and F'{z) = {g{z) - f{z))/{2z) are derived as 



(23) 



(24) 
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